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ABSTRACT 



An unusual Eddington-limited thermonuclear X-ray burst was detected from the accreting neutron star in 2S 0918-549 with the Rossi 
X-ray Timing Explorer. The burst commenced with a brief (40 ms) precursor and maintained near-Eddington fluxes during the initial 
77 s. These characteristics are indicative of a nova-like expulsion of a shell from the neutron star surface. Starting 122 s into the 
burst, the burst shows strong (87 ±1% peak-to-peak amplitude) achromatic fluctuations for 60 s. We speculate that the fluctuations 
are due to Thompson scattering by fully-ionized inhomogeneities in a resettling accretion disk that was disrupted by the effects of 
super-Eddington fluxes. An expanding shell may be the necessary prerequisite for the fluctuations. 
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1. Introduction 

Close to one hundred accreting low-magnetic-field neutron stars 
(NSs) in our galaxy are known to exhibit X-ray bursts from time 
to time 1 (e.g., Cornelisse et al. 2003; Galloway et al. 2008). 
These bursts, lasting roughly 1 min and with a kT « 1 keV ther- 
mal spectrum, are due to thermonuclear shell flashes of hydrogen 
and helium accreted in the top ~ 10 8 g cm -2 (or ~ 1 m) of the 
NS (Woosley & Taam 1976; Maraschi & Cavaliere 1977; Lamb 
& Lamb 1978). The burst recurrence time depends on the rate of 
fuel supply, in other words the mass accretion rate, and the fuel 
mixture (e.g., Fujimoto et al. 1981; Lewin et al. 1993; Bildsten 
1998; Strohmayer & Bildsten 2006). For 96% of all 3387 cases 
studied by Keek et al. (2010), the recurrence time ranges from 1 
hr upward. A few percent of all X-ray bursts are longer than 
usual (e.g., Keek & in 't Zand 2008), due to larger ignition 
depths. These events include "intermediate duration" bursts (in 
't Zand et al. 2005; Cumming et al. 2006; Falanga et al. 2008), 
with durations of order half an hour, and superbursts (Cornelisse 
et al. 2000; Cumming & Bildsten 2001; Strohmayer & Brown 
2002), with durations of order 1 d. These are thought to be due 
to thick (~ 10 10 g crrT 2 ) helium piles in low-M systems and deep 
(~ 10 12 g cm -2 ) ignition of carbon, respectively. 

Helium burning through the 3a process is faster than hy- 
drogen burning. Therefore, if the fuel consists of predominantly 
helium, the energy release is more likely to give rise to super- 
Eddington fluxes. The resulting radiation pressure will expand 
the photosphere which is easily measurable through a combi- 
nation of spectral cooling and increasing emission areas, com- 
bined with a more or less constant bolometric flux. Photospheric 
radius-expansion (PRE) bursts comprise about 20% of all bursts 
(Galloway et al. 2008). Most intermediate-duration bursts, if not 
all, are of the PRE kind. 



For a list, see http://www.sron.nl/~jeanz/bursterlist.html 



The prototype of the intermediate duration burst is from the 
low-mass X-ray binary 2S 0918-549 (in 't Zand et al. 2005). 2S 
0918-549 is thought to be ultracompact in nature with an orbital 
period of 1 hr or less, a hydrogen-deficient donor star to the NS 
and a distance of 5.4 kpc (Nelemans et al. 2004; in 't Zand et al. 
2007). The ultracompact nature was recently supported by the 
tentative measurement of the orbital period of 17.4 min (Zhong 
& Wang 2010). The burst was detected with the BeppoSAX 
Wide Field Cameras and lasted approximately 40 min. Here we 
report a detection of a similar burst from the same system, but 
this time with the much more sensitive Proportional Counter 
Array (PCA) on the Rossi X-ray Timing Explorer (RXTE). We 
identify a new type of X-ray variability in the late stages of the 
burst decay, which we refer to as achromatic late-time variabil- 
ity. We describe the properties of this variability and speculate 
on its origin. 



2. Observations 

Up to May 2010, 2S 0918-549 was observed for 520 ks with 
the PCA (Jahoda et al. 2006), an instrument consisting of 5 co- 
aligned non-imaging proportional counter units (PCUs) which 
combine to an effective area of 6500 cm 2 at 6 keV, have a band- 
pass of 2 to 60 keV with a typical spectral resolution of 20% and 
are collimated to a field of view of 2° x 2° (full-width at zero 
response). No other X-ray sources are known within the field 
of view around 2S 0918-549. Five X-ray bursts were detected, 
two of which have an uncertain thermonuclear origin (Galloway 
et al. 2008). The first four occurred in 2000-2004. On Feb. 8, 
2008, the PCA detected the fifth X-ray burst, when two PCUs 
were active, numbers and 2. The relevant RXTE observation 
identification is 93416-01-05-00. The decay of the fifth burst is 
7 times longer and 20% brighter (in photon count rate per PCU) 
than the next longest/brightest first burst (Jonker et al. 2001). 
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Fig. 1. PCA light curves of long burst from 2S 0918-549 in two band- 
passes and over the 2 active PCUs (top two panels) and the hardness 
ratio between both (third panel). The time resolution is 0.1 s for times 
earlier than 1 s, and 1.0 s for later times. The inset in the first panel 
shows the full-bandpass light curve for the first 60 ms at a time resolu- 
tion of 1 ms with rate in units of 10 4 c s -1 . 

3. Analysis 

3. 1 . Light curves 

Figure 1 presents the light curves of the fifth burst in two band- 
passes and the (hardness) ratio of both. Many aspects of the burst 
are similar to what has been observed before in this and other 
sources. The rise of the burst is very fast. In fact, a precursor can 
be distinguished that merely lasts 0.04 s. This duration is simi- 
lar to that of the record holder (the precursor to a burst from 4U 
0614+09; Kuulkers et al. 2009). The start of the main event is 
1.2 s after that of the precursor. The precursor signifies that su- 
perexpansion happens in this burst which is presumably due to 
an optically thick shell being expelled by strong radiative driv- 
ing (in 't Zand & Weinberg 2010). The time scales are identi- 
cal to the record burst from 4U 0614+09. The hardness peaks 
77 s after burst onset. This identifies the so-called 'touch down' 
point (see spectroscopic analysis below). It coincides with a lo- 
cal minimum (30% below the trend) in the photon count rate for 
energies below 4 keV. Such a decrease is hardly seen above 4 
keV. 

One aspect about this burst is very unusual and is the fo- 
cus of the present paper. At 122 s, a strong variability starts that 
lasts 66 s with ~87% peak-to-peak amplitude. Interestingly, the 
hardness ratio is unaffected: there are no hardness fluctuations 
corresponding to the intensity fluctuations. After this phase, the 
burst decays in an orderly fashion until the end of the observa- 
tion at 313 s when the flux is approximately 10% of the peak 
value. The exponential decay times during this orderly phase are 
149.3 + 3.1 and 1 13.4 + 1.2 s for the low (i.e., below 4 keV) and 
high-energy (above 4 keV) light curves, respectively. 
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Fig. 2. (Top) Full bandpass PCA light curve of a part of the burst that 
exhibits fluctuations, at 1 s resolution. The dashed line shows the ex- 
trapolation of an exponential function fitted to data beyond 190 s. The 
e-folding decay time fitted is 124.8 ± 0.8 s. The inset zooms in on the 
chirp. Two PCUs were active. (Bottom) Fourier power spectrum of same 
light curve, but at 0.5 ms resolution. 



Figure 2 shows a detailed light curve of the variable phase 
and a Fourier power spectrum of the same data. A number of 
features to note here are: 

• the typical time scale is 1 s. There is no variability below 
0.2 s; the Fourier power density above 5 Hz is consistent 
with statistical (white) noise with an upper limit on the frac- 
tional rms of about 1.2% at the frequency resolution applied 
in Fig. 2; 

• A cross correlation analysis between the low (< 4 keV) and 
high energy (> 4 keV) lightcurves shows no time lags in 
excess of 0.01 s; 

• No burst oscillations were detected in any time segment, 
with 3<t upper limits on the fractional rms amplitude (c.f., 
Strohmayer & Bildsten 2006) of 3.9% in 4 s time segments 
and 6.9% in 1 s segments; 

• most dip ingresses are shorter than a few tenths of a second. 
Egresses are almost always longer; 

• dip durations range from 0.2 to 9 s. Dips tend to become 
longer at later times; 

• fluctuations upward are limited to a factor of 1.6 with re- 
spect to the general decay trend (dashed line in Fig. 2). 
Fluctuations downward are down to a factor of 3. The peak- 
to-peak amplitude for the whole period of fluctuations is 
87+1%; 
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Fig. 3. Time-resolved spectroscopy of the burst from 2S 0918-549. The top panel shows the full bandpass light curve of the xenon layers, the 2nd 
panel the bolometric flux (early values suffer from statistical errors between 10 and 20%), the 3rd panel the black body temperature in terms of 
kr eff , the 4th panel the black body radius for a distance of 5.4 kpc, and the bottom panel the reduced^ 2 for the fits per time bin. 



• all upward fluctuations are left and right accompanied by 
downward fluctuations, most clearly for the features at 124 
and 134 s. The converse does not apply; 

• between 138 and 140 s, 4 dips at progressively shorter inter- 
vals are seen, giving the appearance of a chirp. 

3.2. Spectra 

We employed event mode data that have a time resolution of 
125 jus and bin the photon energy in 64 channels. Data from all 
Xenon detector layers and the two active PCUs were combined. 
Spectra were extracted in time bins that are sufficiently small 
to follow the variations in hardness (see Fig. 1), with bin sizes 
measuring 2 ms (until 6 ms after precursor onset; see Fig. 3), 
12 ms (until 1 s), 0.1 s (until 2 s), 1 s (until 200 s) and 5 s (be- 
yond 200 s). The particle and cosmic background were modeled 
per time bin as prescribed by pcabackest with the minimum 
allowed sample time of 16 s. A response matrix was calculated 
through pcarmf version 11.7. Exposure times were corrected for 
detector dead time following RXTE Cook Book 2 . Since the time 
resolution for such a correction is limited to 0.125 s, interpo- 

2 http://heasarc.gsfc.nasa.gov/docs/xte/recipes/pca_deadtime.html 



lations were necessary for the precursor. The spectra were ana- 
lyzed between 3 and 30 keV with XSPEC version 12.5.0 (Arnaud 
1996). A systematic error of 0.5% was assumed per photon en- 
ergy channel. All spectra were corrected for the modeled particle 
and cosmic background. 

A pre-burst spectrum of 2S 0918-549 was acquired from 672 
s of standard-2 data that are available up to 100 s before the burst 
onset. The shape of the pre-burst spectrum is within 3 to 30 keV 
consistent with the broadband 0.1-200 keV spectrum that was 
measured by in 't Zand et al. (2005). The extrapolated 0.1-200 
keV unabsorbed flux is 50% larger at 5.5 x 10~ 10 erg cm _2 s _1 
than the absorbed 3 to 30 keV flux. This is 0.5% of the bolo- 
metric peak flux of the burst (see below). For a canonical 
1.4 M /R=1O km NS, this translates to a mass accretion rate of 
1.0 x 10 16 g s _1 . Allowance was made for the accretion radiation 
in the spectral analysis of the burst data by including a fixed pre- 
burst model, thus indirectly subtracting the background. This 
was not done for the precursor data since it is expected that the 
accretion flow is significantly disturbed (i.e., the flux during su- 
perexpansion is lower than before the precursor). Although the 
accretion radiation was accounted for, it was found that neglect- 
ing this did not make a significant difference for any of the spec- 
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tral results. The accretion flux is always much lower than the 
burst flux for data with a meaningful statistical significance. 

The time-resolved burst spectra were modeled by simple 
black body radiation with interstellar absorption. The absorbing 
column was fixed at JV H = 3 x 10 21 ctrT 2 (in 'tZand et al. 2005), 
which is hardly detectable with the 3-20 keV bandpass of the 
PCA. The results are shown in Fig. 3, in both logarithmic (left 
panel; to focus on the early part of the burst) and linear scales 
(right panel; to focus on the part where the fluctuations occur). 
The bolometric fluence over the measurements of the burst is 
(1.58 + 0.11) x 10~ 5 erg cm~ 2 . Extrapolating the decay beyond 
the end of the observation (through an exponential function with 
an e-folding decay time of 120 s) predicts that 20% of fluence is 
likely to have been missed in the observation. For 5.4 kpc, this 
translates to a total radiative energy output of about 7xl0 40 erg, 
which is consistent with an ignition depth of 4 x 10 9 g ctrT 2 if 
the pre-burst layer is pure helium (see Fig. 8 in in 't Zand et al. 
2005). 

The precursor shows a rapid photospheric expansion driven 
by radiation pressure as the burst emission is super-Eddington. 
Within 30 ms, the photosphere expands to roughly 10 3 km, with 
an average velocity of roughly 0.1c. We do not believe this is 
an optical depth effect, because prior to the burst no gas of suf- 
ficient density is expected up to that height around the NS. It 
would fall to the NS within a dynamical time scale of 1 ms. The 
measured average velocity is roughly 20% of the escape veloc- 
ity from a canonical NS, so this shows that the photosphere is 
strongly driven by radiation for as long as it is visibly expanding. 
Otherwise it would return to the NS already within a height of 
less than 1 km. The photosphere is probably on a shell of thick- 
ness 10 6 " 7 g cm" 2 (c.f., in 't Zand & Weinberg 2010). At the end 
of the precursor the shell radiation moves out of the PCA band- 
pass so that it becomes undetectable. After 1.2 s the burst emis- 
sion reappears. Probably the shell at that time diffused so much 
that it became optically thin and uncovered the underlying radi- 
ating NS. The radiation is still super-Eddington with high tem- 
peratures (>2 keV), moderately increased levels of black body 
radius and a more or less constant bolometric flux. After 70 s 
the flux starts to decrease, signifying the end of the Eddington- 
limited phase. At that time the spectrum shows a hardening last- 
ing approximately 40 s which may be due to Comptonization of 
the thermal photons against a hot receding atmosphere (London 
et al. 1986; Madej et al. 2004). The simultaneous radius mea- 
surements are most probably biased by this systematic effect and 
do not represent true variations of the photospheric size. This 
'touch-down' phase is often seen in Eddington-limited bursts 
(Damen et al. 1990; Galloway et al. 2008; Steiner et al. 2010). 
Subsequently ordinary NS cooling sets in, until the fluctuations 
start at 120 s. As is the case for the hardness ratio in Fig. 1, the 
evolution of the black-body temperature is largely unaffected by 
the fluctuations. The root mean square of the temperature varia- 
tions relative to a fitted exponential function is 1.0 ± 1.6% prior 
to the fluctuations (i.e., between 80 and 120 s) and 2.3 ± 2.3% 
during (i.e., between 120 and 190 s). We conclude that the fluc- 
tuations are achromatic. 



4. Achromatic late-time variability in other bursts 

Achromatic late-time variability was also seen in the previous 
long burst detected from 2S 0918-549, with BeppoSAX-WFC 
(in 't Zand et al. 2005, referred to in § 1). In fact, it looks rather 
similar, but with worse statistical quality (the effective area of 
BeppoSAX-WFC is 45 times smaller than of the PCA). The fluc- 
tuations started 135 s after burst onset, lasted 50 s and the am- 
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Fig. 4. PCA light curve of the part of the superburst from 4U 1820-30 
that shows achromatic variability. For panel explanation, see Fig . 1 . The 
time resolution is 1.0 s throughout. 



plitudes were factors up to 3. Additionally, a very long dip was 
detected starting 200 s after burst onset. The flux decreased by a 
factor of 4 for 70 s. The burst decay is followed during the dip. 
Also for this dip no spectral change could be detected. 

The WFC burst was probably also a superexpansion burst, 
because it shows a rapid decline in black body radius at the onset. 
A definitive identification of superexpansion is not possible since 
a precursor is not detected, but this is not unexpected because 
there would be too few photons in the precursor to be able to 
detect it with the moderately sensitive WFC. The WFC burst has, 
within 20%, the same near-Eddington duration and decay time as 
the PCA burst. These similarities are unlikely to be a coincidence 
and point to the burst as the initiator of the variability. 

Achromatic variability was also detected in two other long 
superexpansion bursts, from SLX 1735-269 (Molkov et al. 2005) 
and 4U 1820-30 (Strohmayer & Brown 2002; Ballantyne & 
Strohmayer 2004; in 't Zand & Weinberg 2010). In the former, 
detected with JEM-X on INTEGRAL, the fluctuations lasted be- 
tween 700 and 1 300 s after burst onset, with a touch down at 400 
s. The burst from 4U 1820-30 was a superburst detected with 
the PCA. Fluctuations in that burst have varying time scales, 
but those at ~ 10 s time scale are also energy independent (see 
Fig. 4). The fluctuations started 3365 s after burst onset and 
lasted until 7200 s, while the touch down for that burst is at 1 100 
s. 

Finally, achromatic variability was seen in a burst from M15 
X-2 (van Paradijs et al. 1990) between 127 and 156 s after burst 
onset, with a touch down time of 94 s. They were attributed to 
dipping activity, although it was not understood why the non- 
burst emission failed to exhibit dips. 

Achromatic variability was not seen in shorter bursts from 
2S 0918-549 that are Eddington-limited but lack superexpansion 
(Jonker et al. 2001 ; in 't Zand et al. 2005). 
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Fluctuations downward were seen in a burst from A 1246-58 
(in 't Zand et al. 2008). A number of intermediate duration bursts 
which (are likely to) have a precursor did not show fluctuations 
(e.g., in 't Zand et al. 2007; in 't Zand & Weinberg 2010). 

In conclusion, based on the aforementioned 6 fluctuating 
bursts, the necessary conditions for fluctuations seem to be the 
presence of superexpansion (i.e., the expulsion of a shell), a burst 
duration in excess of 120 s (although that is susceptible to selec- 
tion effects), and the termination of the Eddington-limited phase. 
Since some long and probably superexpansion bursts lack fluc- 
tuations, another parameter may be at work, perhaps the incli- 
nation angle. Since none of the above bursters show dipping or 
eclipsing behavior in the persistent emission, the inclination an- 
gle must be smaller than i = 70° (White et al. 1995). It seems 
plausible that in bursters which fail to show fluctuations in super- 
expansion bursts, the inclination angle is smaller again, so that a 
disrupted disk does not extend into the line of sight. 

5. Discussion 

We presented an analysis of a type of X-ray intensity variation 
seen in some long intense X-ray bursts with superexpansion. 
Such variations have not been reported previously in such de- 
tail; here we discuss possible origins. 

5. 1 . Origin of fluctuations 

The fluctuations are achromatic and go factors above and below 
the decay trend of the burst. The latter rules out a thermonuclear 
origin, since then the flux would not undercut the decaying trend 
of the initial flash. An example of thermonuclear variability is 
millihertz oscillatory behavior in the non-burst emission of some 
bursters (4U 1608-52 and 4U 1636-536; Revnivtsev et al. 2001; 
Heger et al. 2007; Altamirano et al. 2008). The characteristics 
of these oscillations are different from the achromatic variability 
reported here, both in time scales and spectral behavior. 

The fluctuations downward point to something occasionally 
and partly obscuring the line of sight to the NS. Since the NS 
is small and the dips are « 70% deep, the obscuration appears 
to be due to clouds of optical depth t « 1. The fact that the 
fluctuations are achromatic suggests that the cloud plasma is in 
thermal balance with the radiation field so that the scattering 
of the photons is elastic (i.e., Thompson scattering). The col- 
umn density of the electrons in the clouds, for a uniform den- 
sity, would be N e = t/<x t = 1.7 x 10 24 ctrT 2 (with cr T the 
Thompson cross section). Since the clouds need to be larger than 
the L =20 km size of the NS, the electron density is inferred to 
ben e <NJL* 10 18 ctrT 3 . 

The picture arises that hot plasma clouds move around the 
NS and Thompson scatter the burst photons. If a cloud is in front 
of the NS, burst emission is scattered out of the line of sight. If it 
is at the other side of the NS, but visible to the observer, it is scat- 
tered backward into the line of sight. Since between upward and 
downward fluctuations we see basically unaffected burst emis- 
sion (for instance at 125-129 s, 134-136 s, 162-166 s) only a few 
clouds would be responsible for the fluctuations and there would 
hardly be any sideways scattering. The latter would be consis- 
tent with Thompson scattering, since that is predominantly in 
the forward and backward direction. Still, the flux level between 
fluctuations appears to be somewhat increased. This may be due 
to some sideways scattering or to backward scattering of clouds 
with smaller angular sizes. 

The fluctuations downward are reminiscent of dips as seen 
in the out-of-burst radiation of low-mass X-ray binaries with in- 



clination angles higher than about 70°. Those are thought to be 
due to splashes off the accretion disk that result from the impact 
of the mass stream from the donor star on the disk (Parmar et al. 
1986; Armitage & Livio 1998). The splashes usually populate a 
small range of orbital phases and return to the line of sight pe- 
riodically, obscuring the inner (X-ray radiating) disk, resulting 
in dips in the light curve with typical durations of a few tens 
of seconds. The dip spectrum can be modeled by a combina- 
tion of energy-dependent absorption by cold matter and energy- 
independent absorption by warm photoionized matter with a 
photo-ionization parameter of £ = L/n e r 2 ~ 10 25 ~ 3 erg s~'cm 
(Boirin et al. 2005; Diaz Trigo et al. 2006). The dips reported 
here are not connected to certain orbital phases, but to certain 
phases in long bursts. A simple estimate for the probability that 
the two long bursts occur at the same orbital phase by chance is 
1 min/17.4 min which is only 6%. Also, classical dips in the non- 
bursting emission would also have to occur, but have never been 
reported for 2S 0918-549. Still, our dips are energy-independent 
like part of the classical dips. For log£ S3, L~ 2 x 10 38 erg s -1 
and n e ^ 10 18 cm 4 , the distance of the scattering medium to the 
NS would be of order r % 10 4 km. 

The first down-up-down fluctuation at 124 s appears to be 
the effect of a single cloud that first is in front of the NS, 1 s 
later behind, then again in front and subsequently gone, possibly 
accreted. Almost the same happens at 130 s. This suggests a 2 s 
Keplerian period, which would translate to a typical cloud dis- 
tance from the NS of 3 x 10 3 km which is roughly consistent with 
the photo-ionizing argument above. The chirp at 139 s could be 
explained similarly, with a progressively shorter orbital period. 
Later dips appear to become progressively longer, but peaks stay 
as long and disappear. In the other superexpansion burst detected 
from 2S 09 1 8-549 by the WFC, one dip lasted as long as 70 s. An 
explanation may be that the clouds extend over larger azimuthal 
ranges as time progresses and obscure backscattering. Either the 
clouds become larger in linear size, or they remain as large but 
come closer to the NS. 

Warping of accretion disks may be strengthened by irradi- 
ation from a central source. Previous studies (Petterson 1977; 
Pringle 1996, 1997; Wijers & Pringle 1999; Ogilvie & Dubus 
2001) assumed the central source to be the inner accretion disk, 
but a bursting neutron star would seem viable as well. However, 
the time scales of the radiation are totally different (weeks ver- 
sus minutes) and it is questionable whether burst radiation will 
be as effective in warping the disk. Furthermore, radiative warp- 
ing would occur irrespective of the presence of superexpansion, 
in contrast to what we observe. 

5.2. Accretion disk disruption 

Where would the cloud structure come from? An obvious possi- 
bility is that the clouds arise from material in the inner accretion 
disk. Theoretical investigations have shown that it is possible 
for Eddington-limited X-ray bursts to disrupt the inner accre- 
tion disk. Ballantyne & Everett (2005) treated this issue most 
recently in an attempt to explain spectroscopic measurements 
of the 4U 1820-30 superburst (not the aforementioned fluctua- 
tions); earlier discussions may be found in Fukue (1983), Walker 
& Meszaros (1989), Walker (1992) and Fukue & Umemura 
(1995). Five mechanisms are discussed by Ballantyne & Everett 
(2005), four of which are concerned with the interaction between 
the radiation from the NS with the accretion disk: radiatively and 
thermally driven winds, increased inflows as a result of radiation 
drag (i.e., the Poynting-Robertson effect), and X-ray heating of 
the disk. It is clear that thermally driven winds are ineffective 



6 



in 't Zand, Galloway & Ballantyne: Achromatic late-time variability in thermonuclear X-ray bursts 



due to the large gravitation field near the NS. Radiation pres- 
sure is more effective, but it is not obvious that it can sweep 
the inner accretion because the optical depth of the disk is large 
along the radial travel direction of the burst radiation. Radiation 
drag is probably only important in the innermost region of the 
accretion disk (c.f., Miller & Lamb 1996), at radii smaller than 
20R g « 40 km with R g — GM/c 2 . The final radiative disruption 
mechanism, X-ray heating of the disk surface, may change the 
disk structure through for example puffing up the surface on a 
dynamical time scale and settling it on viscous time scales. A 
non-radiative fifth mechanism is concerned with the mechanical 
interaction between the radiatively driven outflow and the ac- 
cretion disk. This effect was dismissed by Ballantyne & Everett 
(2005) for a radiation-driven wind on the basis of that carrying 
insufficient kinetic energy to overcome the potential energy of 
the wind. However, the circumstances may be different for su- 
perexpansion bursts where the shell carries much more density 
than the previously considered disk. In §5.2.1 and 5.2.2 we con- 
centrate on the two most likely mechanisms for disrupting the 
disk. 

5.2.1. Disruption by shell 

The geometrically thin but optically thick shell expands at a ve- 
locity of order 10% of the speed of light. The column thick- 
ness of the shell at launch is at maximum ~1% of the igni- 
tion depth or y 10 6 ~ 7 g crrT 2 and diffuses as R 2 . The ra- 
dial profile of the density of the accretion disk at mid plane, 
for a gas-supported disk (as applicable for an accretion rate of 
~1% of Eddington) with electron-scattering opacity (Shakura 
& Sunyaev 1973; Svensson & Zdziarski 1994; Ballantyne & 
Everett 2005) and assuming a Gaussian vertical density profile, 
is given by 

p = po (R/R g y 3V2() J(R) 2/5 g cm" 3 (1) 
with R g = 2M\a km (M\ A is the NS mass in units of 1.4 M Q ), 

J(R) = 1 - ^6R g /R (2) 
and 

po = 47 ^M 16 2 W'» g cm- 3 (3) 

(ds.A is the distance in units of 5.4 kpc, M\$ the accretion rate in 
units of 10 16 g s _1 and oro.i the disk viscosity parameter in units 
of 0.1). The mean molecular weight is assumed to be /i = 1.3. 
For a = 0.1, Mi 6 between 1 and 10 (for radiation efficiencies of 
i] = 1 and 0.1, respectively) and y = 10 6 g crrr 2 , the mass col- 
umn swept up by the diffusing shell will have reached the shell's 
own diminished column mass at 50 to 150 km above the NS sur- 
face. For y = 10 7 g crrr 2 this range becomes 1 10-440 km. We 
assume here that the disk is magnetically truncated at 25M~ 2 ^ 7 
km from the NS (for a canonical NS with B = 10 8 G; this does 
not have to be accompanied by pulsations, see Frank et al. 2002). 
Away from the mid plane the disk density is lower. At 3-4 scale 
heights above or below the disk (this is approximately equiva- 
lent to the t = 1 level), the shell may be able to disrupt the disk 
up to a distance of ~ 3000 km. Thus, the shell appears able to 
blast away a small radial part of the disk, out to perhaps a few 
tens of NS radii, and ablate the surface of the disk up to dis- 
tances of IO 3 - 4 km. Note that with a suggested orbital period of 
17.4 min (Zhong & Wang 2010), the accretion disk has a size of 
1.7xl0 5 km. 



It is difficult to determine whether the energy and momen- 
tum are conserved in the the disk/shell interaction. The appar- 
ent expansion velocity of the shell is only ~20% of the escape 
velocity at the NS surface which allows to bridge only 4% of 
the potential energy required to escape, while the shell is ob- 
served to bridge 99% (i.e., the maximum photospheric radius is 
10 3 km). Obviously, the transformation of radiation to kinetic 
energy is a dominant term in the energy balance. An assess- 
ment of that is beyond the scope of the present observational 
paper. Conservation of momentum is roughly upheld. The ram 
pressure of the shell at 100 km is pv 2 = y(R/R NS ) 2 (0.lc) 2 /D a 
10 17 dyne cm 2 with D = 1 m the thickness of the shell or per- 
haps a few orders of magnitude smaller if the shell thickness 
grows with the shell radius. The density of the disk at 100 km 
is of order 10~ 3 g crrr 3 (Eq. 1). From M = 2nRHpv this results 
in a radial inward velocity of « 20 km s _1 . The ram pressure of 
the disk 3 then is « 10 7 dyne cm 2 at 100 km. However, the gas 
pressure of the disk is larger and will take over counteracting the 
shell. For an ideal gas it is ~ 10 15 dyne crrr 2 for a disk temper- 
ature of 10 7 K. This order-of-magnitude calculation shows that 
the shell may have enough momentum to sweep the inner 100 
km of the accretion disk. 

5.2.2. Disruption by radiation pressure 

In the previous section it is shown that in the first 10 _1 s the 
shell may disrupt the accretion disk over the inner 10 3 ~ 4 km, 
possibly even clearing it completely over the inner 10 2 km. 
Thus, an environment is created which is more susceptible to the 
near-Eddington fluxes that ensue for 70 s. Ballantyne & Everett 
(2005) modeled the effects of super-Eddington fluxes on the ac- 
cretion disk without the preceding shell disruption, but with oth- 
erwise similar circumstances (same burst flux that stays near- 
Eddington for at least 70 s; same type of binary orbit). They find 
that the gas in the surface layers of the disk with n e = 10 17 cm -3 
is accelerated radially to «0.1 c within at most » 30 s. Thus, the 
radiation is likely to further clear the parts of the disk already 
disrupted by the shell or, at least, prevent it from falling back. 
The effect may even be larger, because the shell will likely puff 
up the disk, increasing the solid angle of accretion disk material 
as seen from the NS so that a larger portion of the burst luminos- 
ity will be involved in the radiative disruption. 

5.3. Accretion disk resettlement 

If the disk is disrupted, it is plausible that it deviates considerably 
from cylindrical symmetry. The push by the shell and its lateral 
expansion will likely puff up the disk, as will the radiation pres- 
sure due to turbulence. If the effect of super-Eddington radiation 
subsides, the accretion disk is expected to return to the immedi- 
ate surroundings of the NS, because the feeding of the disk by 
the donor star continues. We speculate that it is this return of the 
disk that is responsible for the achromatic fluctuations. 

How can, then, the delay time be explained of 45 s between 
the end of the Eddington-limited phase and the onset of achro- 
matic variability? The free fall time over 10 3 km of 0.5 s is much 
shorter, so it involves something else than a straightforward fall 
back of disk material. Could it be the viscous dissipation time 
scale of the disk? This time scale is, for a gas-pressure dominated 
accretion disk, given by (Svensson & Zdziarski 1994; Ballantyne 

3 It is noted that the larger disk ram pressures calculated in Frank 
et al. (2002) are for azimuthal motion instead of radial, as is applicable 
for balancing magnetic pressure 
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& Everett 2005) 

fvisc = 4a Q fM^M^ /5 (R/R g ) 7/5 J(Rr 2/5 s. (4) 

At 100 km, f v ; sc ranges between 60 and 960 s. This is in line with 
the measured delay time of 45 s. However, not so for the 2000 s 
delay time for the superburst from 4U 1820-30. Mi 6 is 11 times 
larger and f v j sc 2.5 times smaller while the measured value is 50 
times larger. 

Rather than the time it takes the disk to respond to changes 
in the environment, the onset of the fluctuations may be deter- 
mined by the level of irradiation (see § 5.2.2). The radiation 
pressure is equal to L/4-7tR 2 c. The disk will be disrupted as long 
as L > AnR 2 cP mm with P mm the (uncertain) ram pressure of 
the in-flowing accretion disk. For L — 2 x I0 3& d^ 4 erg s _1 and 
R = 100 km, P rdm would have to be 5 x 10 13 dyne ctrT 2 . The 
uncertainties are large in the disk parameters (M, T, R, P nm ) 
and the physics employed may be oversimplified. Therefore, it 
is difficult to judge the validity of this effect. 

The duration of the achromatic variability, 66 s, may be re- 
lated to the radial extent over which the clouds reside when they 
start to fall back. The longer the near-Eddington phase lasts, the 
further away disk material is blown off by the radiation pressure. 
Indeed, there is a relation between the duration of the variability 
and that of near-Eddington fluxes of the cases discussed in § 4 
(roughly linear). 

6. Summary and outlook 

We report achromatic late-time variability in thermonuclear X- 
ray bursts from 2S 0918-549 and a few other low-mass X-ray 
binaries, and suggest it may be explained by the settlement of 
the accretion disk near the NS that is disrupted by the combi- 
nation of a shell expulsion and a near-Eddington flux. Order- 
of-magnitude calculations of the effect of the various disruption 
mechanisms suggest that this explanation is plausible, although 
there is a need for follow up by numerical time-dependent cal- 
culations of the dynamical response of the accretion disk to the 
expanding shell and wind, the near-Eddington radiation and the 
subsequent settling disk. We note that such a numerical inves- 
tigation may also be interesting for the study of less energetic 
classical novae on massive white dwarfs where the shell is ex- 
pected to be less massive so that the disk may survive the explo- 
sion (Sokoloski et al. 2008). Besides these theoretical investiga- 
tions, detecting more bursts with this variability with a variety of 
burst durations would be constraining, particularly from X-ray 
binaries with known orbital periods. An observation that would 
reject our hypothesis would be achromatic variability in bursts 
without superexpansion with a device that is sensitive enough to 
detect 10 ms precursors. Finally, we note that our findings imply 
a warning: not all features in X-ray burst light curves originate 
on the neutron star (see also Shaposhnikov et al. 2003) and in 
some cases these additional effects can be significant. 
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